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SUMMARY 


This report is presented in fulfillment of the statement of work of contract 
No. NASI -6969 for Phase II as amended on May 13, 1968. The statement of work 
required that TRW analyze the effects of Scout error sources, both individually 
and collectively, on the final trajectory conditions of the following two missions: 

(1) Reentry mission. 

( 2 ) Five stage escape mission with both fourth and fifth stages spin stabilized. 


TRW was to use the same methods and error sources developed under Phase I of the 
subject contract, with the Government supplying fifth stage data and error sources 
for the escape mission and new pitch profiles for both flights. 


The results of analyzing the reentry mission produced the following one sigma 
errors at the time of nominal 400,000 feet altitude : 


nominal 

conditions 

Velocity 26202 ft/ sec 

Flight Path Angle -3*06 deg 

Inclination 50.425 deg 


mean 

dispersion 

- 15.29 ft/sec 
.0736 deg 
.047 deg 


standard 

deviation 

87.0 ft/sec 
•553 deg 
.432 deg 


The analysis of the five stage escape mission resulted in the following dis- 
persions at fifth stage burnout : . 


Velocity 

Altitude 

Flight Path Angle 
Inclination 


nominal 

conditions 

39199.7 ft/sec 
IO 7.56 nmi 
11.769 deg 
37.667 deg 


mean 

dispersion 

2.02 ft/sec 
.424 nmi 
-.007 deg 
.0168 deg 


standard 

deviation 

78.7 ft/sec 
5 .O 3 nmi 
•377 deg 
.0455 deg 


INTRODUCTION 

This report is presented in fulfillment of Phase II of contract No. NASI -6969 
as amended on May 13 , 1968. The statement of work called for the analysis of two 
missions : 

(l) A four stage trajectory designed to reenter the Earth's atmosphere. 

TRW was instructed to determine the expected dispersions in the tra- 
jectory at an altitude of 400,000 feet. 
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(2) A five stage mission with the fifth stage spin stabilized designed to 
enter interplanetary space. The data for the fifth stage was supplied 
by the Scout Program Office at Langley Research Center. 

The error source magnitudes and the method of analysis were to be the same as 
developed and utilized in Phase I of the subject contract. A brief description of 
these methods follows but the interested reader is referred to the Phase I final 
report (Reference l) for a more detailed description of the scout vehicle. 

The Scout simulation equations were assumed to be a quadratic form using 
coefficients previously calculated in the TRW N-Stage Trajectory Program (MVNS). 
There were three component parts: 


Linear Errors : 





linear 


( 1 ) 


where X^ 


% 


is the state vector at a stage time 

is the vector of linear sensitivity coefficients associated with 
.. .th 

the 1 error source 


is a random number with standard deviation equal to c\ 


Non-linear Errors: 



( 2 ) 


where X^ 


is the non-linear state vector 


VV 


is the piece-wise linear function for the non-linear source j , 
formed by its effect at ± la, ±2a, ±3a. 


Cross-term Errors: 


^CT 



12 

I 

t=k+l 


% A k \ 


( 3 ) 


where X^ 


Ski 


is the cross-term vector 
is the cross-term sensitivity 


finally : 


X 


-L + ^NL + ^CT 


(*) 
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The procedure is divided into three consecutive steps: (l) the generation 

of sensitivity coefficients of the state vector at the stage ignition times and 
at burnout or reentry with respect to error sources deviations; (2) the generation 
of cumulative distribution functions of desired orbital and reentry parameters 
using the sensitivity coefficients and the statistic of the error sources; ( 3 ) the 
analysis of these distribution functions. In Figure 1, a block diagram describes 
the relationship between the TRW programs which was used in the analysis. 

Each step will be discussed separately in some detail. 


Sensitivity Coefficient Generation 


The state of any system, as a function of time, may be considered as a 
function of initial state as well as a number of performance parameters (thrust, 
commanded turning rates, aerodynamic, properties , etc.). The effects caused 
by nominal values of these parameters serve to define a nominal state at some 
time of interest, t k , i.e. , 


= f (X o' P ' *0* V 


( 5 ) 


However, if these parameters are allowed to depart from their nominal values, 
a perturbed state vector is generated: 


X (t k ) = f (X , P + SP, t Q ; t k ) 


( 6 ) 


If the difference in the two state vectors is expanded into a Taylor series 
about the nominal value: 

-3s , . 

6X i - L % 6P j * 2 L C «k 6P j » k * - 1 " X - 2 > ••• 6 (T) 
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where : 


6x i = 


the component of the difference between perturbed and 
nominal state vectors 


C -k 

ajk 


is the dimension of the performance variation vector 
is the partial derivative dX^/dP^ 

is the second partial derivative 9 X^/dP^F^ 


The number of terms in the expansion is sufficient to adequately represent 
the functional dependence of the state variation with respect to these parameters. 
If this expansion is truncated to include only second order terms, there are 
q first order parameters and second order coefficients. This would 

necessitate , c * computer runs which for the values of q modelled for the 

Scout analysis (~100) is unrealistic. 


However, it is certainly true that certain error sources will contribute 
more to the state vector dispersion than others and only these major contributors 
need be investigated for nonlinear effects. Thus, if the linear effect of an 
error source on the dependent variables was found to be negligible in relation 
to that of more significant error sources, no attempt was made to investigate 
such minor sources in any greater detail. This procedure has been found by 
TRW Systems to be most efficient, since the time spent in modeling any error 
sources is roughly proportional to its overall contribution to the output. The 
procedure for generation of sensitivity coefficients was as follows: 

1. Using the MVNS simulation of the nominal mission, perturb the 
system with a +3 sigma variation for each error source. Take 
the square root of the sum of the squares of the deviations and 
compare each deviation with the RSS total. If the individual 
deviation is less than 0.3 of the RSS total, the effect of the 
error source will be represented by a linear sensitivity function. 

This criterion guarantees that the error source deviation, when 
squared has a contribution of magnitude less than the total sum 
of the squares. 
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2 . 


The selected sources were perturbed by ±lcx and ±2o and 4-3 o to 

determine their degree of nonlinearity. Both postive and nega- 
tive perturbations are required, since it is possible for the curve 
to be symmetrical while still nonlinear. If it were found that the 
sensitivity curve had a nonlinear shape, a straightline approximation 
was used to represent the partial as shown in Figure 2. 

3 . The next relationship to be found was the cross-correlations 
between these most significant variables. The difficulty of 
this problem is increased by the fact that the important combina- 
tions are not all known a priori. TRW Systems isolated and 
identified several combinations for the Phase I reference mission, 
so that it was not necessary to use the corresponding large amounts 
of machine time for each subsequent mission. 

The non-linear sources investigated in this analysis were thrust misalign- 
ments and coning rate errors, since not only were these sources large contributors 
to the combined errors, but it is clear that major effects such as velocity loss 
will be the same whether the misalignment is positive or negative. Conversely, 
a large effect such as specific impulse tends to be linear, more impulse 
produces more velocity and vice-versa. In addition, this type of error source 
has no out of plane or in plane directional effects. To summarize, the non- 
linear sources used in the analysis of the reentry mission were: 

Thrust misalignments, first stage (TMPl) 

Thrust misalignments, first stage (TMfl) 

Thrust misalignments, second stage (TMP2) 

Thrust misalignments, second stage (TMI2) 

Thrust misalignments, third stage (TMP3) 

Thrust: misalignments , third stage (TMI3) 

Coning impulse, fourth stage (WC4P) 

Coning impulse, fourth stage (WC4T) 

In the escape^ mission the following two sources were added to the list: 
coning impulse, fifth stage (WC5P), coning impulse, fifth stage (WC5Y) 

The cross correlations used in the Scout analysis were those combining thrust 
misalignments at second and third stages with their respective control system 
errors. These combinations are necessary to discover any of the effects 
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of control system errors since if they are operating on a trajectory with no 
upsetting moments as second and third stage are, there will be no reason for 
any control system activity and the nominal burnout values will be reproduced. 
Therefore it was necessary to provide a large and relatively constant upsetting 
force which led logically to the choice of thrust misalignments. The pairings 
which were used in "both missions were: 


Second Stage: 


Thrust Misalignment, Pitch! 
Thrust Misalignment, Pitch: 
Thrust Misalignment, Pitch: 
Thrust Misalignment , law: 
Thrust Misalignment, Yaw: 
Thrust Misalignment, Yaw: 

Third Stage: 

Thrust Misalignment , Pitch: 
Thrust Misalignment, Pitch! 
Thrust Misalignment, Pitch: 
Thrust Misalignment, Yaw: 
Thrust Misalignment, Yaw: 
Thrust Misalignment, Yaw: 

It can be assumed that an output 
variables by the relation 


Rate Gain, Pitch 
Dead Band, Pitch 
Roll Offset 
Rate Gain, Yaw 
Dead Band, Yaw 
Roll Offset 

Rate Gain, Pitch 
Dead Band, Pitch 
Roll Offset 
Rate Gain, Yaw 
Dead Band, Yaw 
Roll Offset 

variable is related to two independent 


5w = f (5x,Sy) = C x &x + C y By + C xx &x 2 + + C^bxSy 

&w = C x &x + C xx &x 2 + Cy&y + CyyBy 2 + C^BxSy 

sy * S 

f x (x) f 2 (y) 

6w - f^Cx) - f 2 (y) 

C xy “ &x8y 


( 8 ) 

(9) 


( 10 ) 
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In this equation f-^ (x) and f 2 (y) are represented by the piecewise linear func- 
tions described in the previous paragraph. 

The effect of carrying out the three-step procedure outlined above has 
to reduce the nonlinear process (i.e. , the system and its environment) to a 
polynomial approximation . This procedure reduced machine running time per 
sample from 2 to 3 minutes on the TRW Systems N-Stage program to 2 seconds 
on the Monte Carlo Error Analysis Program, while still including all significant 
non-linearities . 

This polynomial fit was used to relate the input quantities to the state 
vector at burnout , and any unpowered propagations to other points on the 
trajectory were calculated by an analytical ephemeris generator. 

Monte Carlo Analysis 

Using the nomenclature of the previous section, let us assume that there 
are q error sources, of which p are found to have either nonlinear, cross- 
coupling terms or both. Since the nonlinear terms have +1 CT , +2 a and +3cr 
deviations, each nonlinear error source will have a 6 x 6 matrix of sensitivity 
coefficients, i.e., the value of the output state vector as a function of +n cr . 

These values will be used to construct a piece-wise linear function as seen 

( 2 ) 2 

in Figure 3 and replace the square terms C: ' 6 P. in Eq. 7 . The advantage of 

1J-K J 

this approach is that it is not necessary to fit an arbitrarily truncated 
power series when the curve may contain higher components . 

For each cross coupling pair of error sources, there will be a single 
coefficient for each output state vector component or a 6x1 vector for the 
output state vector. The error sources which have only. linear terms can all 
be compressed into a single equivalent variable by the following steps. Given 
that ; 

X = [C] P (11) 
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where 


X is the output 6x1 state vector, 

P is the linear error source vector of dimension q-p , 

C is the 6 x (q-p) matrix of sensitivity coefficients, 

T T T 

then, XX = CPF C . 

Taking the expectation of each side gives 

S x = CSpC T (12) 

where 

E x is the covariance matrix of the output state vector due to linear error sources; 

£p is the covariance matrix of the linear error sources# 

Thus, all the linear error sources and their statistics can be compressed into a 
single 6x6 covariance matrix. 

The Monte Carlo analysis progresses through the following steps: 

(1) Using the equivalent linear covariance matrix and a random vector 
generator, choose a random output vector from this covariance matrix £ x • 

(2) Using the statistics of each nonlinear source, choose a random value of 
this error source. Enter the piecewise linear function associated with 
the source and calculate the value of the output state vector dispersion 
due to this nonlinear effect. 

(3) Using the statistics of each pair of cross coupling error sources 
choose pair random, values. Multiply the product of these values 
by the cross coupling sensitivity coefficient vector to calculate the 
dispersion due to the combined effects of these sources. 

(4) Add the three dispersions due to linear, nonlinear and cross coupling 
effects to the nominal state vector at burnout to calculate one sample 
of the analysis. Using this perturbed vector calculate the variables 
of interest. For the Scout analysis, the following variables were 
calculated and compared with their nominal values: 

a. Geocentric Radius Vector 

b. Inertial Velocity 

. : ''V» 
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c. Inertial Flight Path Angle 

d. Relative Velocity 

e. Flight Path Angle w/r to the Air Mass 

f . Semi-major Axis 

g. Eccentricity 

h. Inclination 

i. Longitude of the Ascending Node 

j . Argument of Perigee 

k. Downrange Distance from Launch Site 

l. Apogee Altitude 

m. Perigee Altitude 

n. Period 

o . Longitude 

p . Latitude 

q . Altitude 


Statistical Analysis of Samples 

The analysis of the sample statistics is relatively straightforward. 

The samples of any variable are used to calculate a mean and second central 
moment. The samples are then ordered by size and a cumulative distribution 
Is plotted. This curve shows the correct percentile levels for the variable 
•>ven though the distribution may be highly non-Gaussian. 

The equations which are used in the Statistical Processor program (PROC) 
ill be presented here. The same equations apply to both scalars and vectors 
unless a distinction is made in the description of the sample statistic. All 
computations are performed in double precision. 

AVERAGE OR MEAN VALUE 

The sample mean of a random vector X, denoted by X, is computed by PROC 
using the following equation: 


N 



i=l 
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In this and all following equations, N is the total number of samples used 

til 

in the computation and X^ is the value of the random vector X for the i 
sample of the Monte Carlo simulation. 

COVARIANCE MATRIX AND CORRELATION MATRIX 

The covariance matrix of a random vector X, denoted by EL-, is defined 
by the following equation in which E indicates the expectation operator. 

E 7JC = E(X - X)(X - X) T 

This equation can be expanded and rewritten as follows: 

E^ = E(XX T ) - (X)(X) T 


The equation used by PROG to compute a sample covariance matrix is similar 
and can be written 


2 XX 


1 

N-l 


N 


f N \ I 

N f 

s x.x! 

i=l 1 1 

J. 

N 

S X 1 
\i»l j 1 

S x 

i ljsl / J 


CUMULATIVE DISTRIBUTION FUNCTION 

A cumulative distribution function (CDF) is the probability that a random 
variable, X, is less than a specified number. This information is displayed 
by PROC as a graph of the probability that the random variable is less than the 
specified number versus the specified number. This is accomplished by ordering 
all of the samples in the order of increasing numerical value and plotting the 
percentage of the samples that are less than a given value, say v, as a function 
of v . 


12 



In Table 1 are listed all the error sources and their magnitudes as they were 
Listed in Reference 1. The items in part I of the table were not the subject of 
any investigation in Phase I of the contract. Those in part II were the results of 
an investigation of log book and test data performed by TRW and described in 
detail in Reference 1. The fifth stage error used in the escape mission were 
supplied by the Scout Office in a letter. 
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TABLE 1. ERROR SOURCE LIST 


I. Errors From Previous Investigations 

Description 

IA First Stage Errors 

1. Propellant Weight Uncertainty 

2. First Stage Inert Weight 

3. Second Stage Inert Weight 

4. Third Stage Inert Weight 

5. Fourth Stage Inert Weight 

6. Specific Impulse 
7- Mass Flow Rate 

8. Thrust Misalignment - Pitch 

9, Thurst Misalignment - law 

10. Drag Coefficient (CAO) 

11. Normal Force Coefficients (CNa) 

12. Normal Force Coefficients (CN6q) 

13. Side Force Coefficients (CyP) 

14- Side Force Coefficients (Cy6r) 

15. Roll Moment Coefficients (CL6p) 

16. Roll Moment Coefficients (dtp) 

17. Pitch Moment Coefficients (Cmo 


18. 

Pitch Moment Damping Coefficients (Cm6q) 

CMDQ 

.002 

19. 

Pitch Moment Coefficients (Cmq) 

CMQl 

.002 

20. 

Pitch Moment Coefficients (Cma) 

CMAL 

.05 

21. 

Yaw Moment Coefficients (Cnp) 

CNBA 

»33 

22. 

Yaw Moment Coefficients (Cn5r) 

CNDR 

.33 

23. 

Yaw Moment Damping Coefficients (Cnr) 

CNR1 

.33 

24 . 

Density Variation 

DRHO 

*0667 

25. 

Wind Profile 

FWN1 

- 

26. 

Jet Vane Drag Coefficient 

CDV1 

.1 

27. 

Jet Vane Drag Coefficient 

CDV2 

.1 

28. 

Jet Vane Drag Coefficient 

CDV3 

.1 

29. 

Jet Vane Side Force Coefficient 

LDA2 

.1 


Where no units appear, the magnitude is 

a fraction of nominal. 



** This error source magnitude was inadvertently input to Both analyses 
as .002. The effect of the error source was so small, however, that 
it could not effect the results even if the la value were increased 
25 times. 


la Mag. 


Symbol 

Mag . 

PW10 

.0006 

SIW1 

.0083 

SIW2 

.0041 

SIh[3 

.00093 

SIW3 

.0024 

ISP1 

.0018 

MFR1 

.014 

TMP1 

1.67 mrad 

TMY1 

I .67 mrad 

CA01 

.01 

CNAL 

.2 

CNDQ 

.2 

CYBA 

.33 

CYDR 

.33 

CLDP 

.1 

CLP1 

.1 

CM01 

.05** 



la Mag.. 


Description 

30. Roll Moment due to Yaw Axis Shift of Static 

Symbol 


Margin 

31. Roll Moment due to Pitch Axis Shift 

of Static 

LSMY 

.01 

Margin 

32. Pitch Moment due to Roll Axis Shift 

of Static 

LSMP 

.01 

Margin 

33* Pitch Moment due to Yaw Axis Shift 

of Static 

MSMR 

.1 

Margin 

34. Yaw Moment due to Pitch Axis Shift 

of Static 

MSMY 

.1 

Margin 

35. Yaw Moment due to Roll Axis Shift of Static 

NSMP 

.1 

Margin 


NSMR 

.1 

IB Second Stage Errors 

36. Propellant Wt Uncertainty 


PW20 

.00054 

37* Specific Impulse 


ISP2 

.00094 

38. Mass Plow Rate 


MFR2 

.01 

39* Thrust Misalignment - Pitch 


TMP2 

1.67 mrad 

40. Thrust Misalignment - Yaw 


TMY2 

1.67 mrad 

41. Second Stage Aerodynamics (CDO) 


CD02 

.1 

42. Second Stage Aerodynamics (Cm.) 


CM2 

.1 

43* Second Stage Aerodynamics (§) 


ZET2 

.1 

IC Third Stage Errors 

44- Propellant 'Weight Uncertainty 


PW30 

.0006 

45 * Specific Impulse 


ISP3 

.0014 

46. Mass Flow Rate 


MFR3 

.018 

47* Thrust Misalignment - Pitch 


TMP3 

.557 mrad 

48. Thrust Misalignment - Yaw 


TMY3 

.557 mrad 

ID Fourth Stage Errors 

49* Propellant Weight Uncertainty 


PW40 

.00034 

50. Specific Impulse 


ISP4 

.006 
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Description 

Symbol 

Irr Mag. 
Mae. 

51. 

Mass Flow Rate 

MFR4 

.018 

52, 

Thrust Misalignment - Pitch 

TMP1 

.5 mrad 

53- 

Thrust Misalignment - Yaw 

TMY1 

•5 mrad 

54. 

Coning Rate - Pitch 

W4CP 

.03 rad/sec 

55. 

Coning Rate - Yaw 

W4CY 

.03 rad/sec 


II'. Error Sources Resulting from Investigation 


IIA First Stage Errors 

Symbol 

la Mag. 

56. 

Proportional Gain Error - Pitch 

KPP1 

.0243 

57. 

Proportional Gain Error - Yaw 

KPY1 

.0306 

58. 

Proportional Gain Error - Roll 

KPR1 

,0233 

59. 

Rate Gain Errors - Pitclj 

KRP1 

.044 

60. 

Rate Gain Errors - Yaw 

KRY1 

.021 

61. 

Rate Gain Errors - Roll 

KRR1 

.0262 

* 62. 

Random Uncompensated Gyro Drift - Pitch 

DTEP 

.125x10 ^ rad/sec 

63. 

Random Uncompensated Gyro Drift - Yaw 

DTEY 

. 200X10 - ^ rad/sec 

64. 

Random Uncompensated Gyro Drift - Roll 

DTER 

.328X10 ^ rad/sec 

65. 

Mass Unbalance Roll Gyro - Input Axis 
/ 

KRIA 

1.03 xlO ^ rad/sec/ 
ft/sec 2 

66. 

Mass Unbalance Pitch Gyro - Spin Axis 

KPSA 

1.08 XlO - ”* 7 rad/sec/ 
ft/sec 2 

67. 

Mass Unbalance Yaw Gyro - Input Axis 

KYIA 

1.03 xlO ^ rad/sec/ 
ft/ sec 2 

68. 

Anisoelastic Errors - Pitch 

KPAN 

3 .13xl0~' 1 ^' rad/sec/ 
ft^sec^ 

69. 

Anisoelastic Errors - Roll 

KYAN 

3. 13 XlO -11 rad/sec/ 
ft 2 / sec^ 

70. 

Anisoelastic Errors - Yaw 

KRAN 

3. 13 xlO -11 rad/sec/ 
ft 2 / sec* 
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Symbol 

la Mag. 

71. 

Rate Gyro Bias - Pitch 

DPBE 

3-57 mrad/sec 

72, 

Rate Gyro Bias - Yaw 

DYBE 

3.57 mrad/sec 

73. 

Rate Gyro Bias - Roll 

DRBE 

3.57 mrad/sec 

74. 

Rate Gyro Misalignment - Yaw 

TIRG 

' 1.45 mrad 

75. 

Rate Gyro Misalignment - Roll 

TRRG 

1.45 mrad 

76. 

Vehicle Alignment Errors - Pitch 

THOP 

5.76x10 rad 

77. 

Vehicle Alignment Errors - Yaw 

THOY 

5.43x10 rad 

78. 

Vehicle Alignment Errors - Roll 

THOR 

5.25xl0“ 5 rad 

79. 

Intervalometer and Torque r Scale Factor ■ 

DKSG 

.0035 

Fin 

Misalignment Errors 



80. 

Normal Force Error/Yaw Fins 

CNDR 

.576x10“^ rad 

81. 

Side Force Error/Pitch Fins 

CYDQ 

.576x 10" 4 rad 

82. 

Roll Moment Error/Yaw Fins 

CLDR 

.576x10”^ rad 

83. 

Roll Moment Error/Pitch Find 

CLDQ 

.576x10“^ rad 

84* 

Pitch Moment Error/Yaw Fins 

CMDR 

. 576 x 10 "^ rad 

85- 

Yaw Moment Error/Pitch Fins 

NCDQ 

•576x10”^ rad 

86. 

Timer Error - First Step 

TIML 

.078 sec. 

87. 

Timer Error - Second Step 

TIM2 

.004 sec. 

88. 

Timer Error - Third Step 

TIM3 

.003 sec. 

89. 

Timer Error - Fourth Step 

TIM4 

.003 sec. 

IIB. Second Stage Errors 



90. 

Rate Gain Error - Pitch 

KRP2 

.044 

91. 

Rate Gain Error - Yaw 

KRY2 

.026 

92. 

Dead Band Error - Pitch 

DBP2 

.1 

93. 

Dead Band Error - Yaw 

DBY2 

.1 

94. 

Roll Offset 

R0E2 

.25 deg. 

95. 

Control Motor Misalignment - Pitch axis 

\ 

C2PY 

,0033 deg. 

96. 

Control Motor Misalignment - Yaw axis ! 

C2IP 

.0033 deg. 

97. 

Timer Error - Sixth Step 

TIM6 

.003 sec. 

98. 

Timer Error - Seventh Step 

TIM7 

.003 sec. 


17 




99. 

Rate Gain Error 

- Pitch 

i KRP3 

.044 

100. 

Rate Gain Error 

- law 

KRY3 

.021 

101. 

Dead. Zone Error 

- Pitch 

DBP3 

.1 

102. 

Dead Zone Error 

- law ; 

DBY3 

.1 

103. 

Roll Offset 


R0E3 

.25 deg. 


IID. 

Fifth Stage Errors 

Symbol 

la Mag* 

104. 

Mass Flow Rate 

MFR5 

.018 

105 . 

Propellant Weight Uncertainty 

PW50 

.00033 

106. 

Specific Impulse 

ISP5 

.0015 

107 . 

Fif^h Stage Inert Weight 

SIW5 

.00033 

108 . 

Thrust Misalignment, Pitch 

TMP5 

.667 mrad 

109 . 

Thrust Misalignment, Yaw 

TMY5 

.667 mrad 

110. 

Coning Rate, Pitch 

W5CP 

.0463 rad/sec 

111. 

Coning Rate, Yaw 

WSCY 

.0463 rad/sec 


18 



REENTRY MISSION ERROR ANALYSIS RESULTS 


The mission input data is received from the Scout Program Office and the 
nominal trajectory which resulted is presented in Table 2. 

Table 3 lists the mean, variances and extreme values of the output parameters 
and Table 4 is a summary of the most significant error sources and their three 
sigma contributions to the dispersions in range velocity, flight path angle and 
inclination. 

Tables 5, 6, 7 and 8 are the individual three sigma contributions of each 
error source to altitude dispersion, velocity, flight path angle and inclination 
measured at each stage ignition and the nominal 400,000 ft time. 

Tables 9 and 10 present the non-linear and cross-term data which went into 
the Monte Carlo results. In Table 9, there is presented the non-linear data for 
eight different sources. The position and velocity deviations in a radial, down- 
range and out-of-plane coordinate system were calculated for an error source input 
of +3c, +2a, +la, -la, -2o and -3a. The Monte Carlo program then interpolated 
between these points. 

Table 10 displays the cross-term data involving the effects of a pair of 
error sources acting simultaneously with both sources at their 3a level. The sources 
chosen were all in the control system which will not respond except in the presence 
of a disturbing force. Thrust misalignment, being the largest disturbing force 
was chosen to excite these errors. The deviations are represented in the same 
radial, downrange and out-of-plane coordinate system and they have been normalized 
by dividing by nine times the product of the two standard deviations. Analysis of 
these results by comparing the deviations with the 3a thrust misalignment deviations 
indicates no significant cross term effects. 

Figures 3-19 depict the cumulative distribution function of each of seventeen 
output variables. 

The results indicate that the velocity dispersion is approximately 15 ft/sec 
greater than in the polar case examined in Phase I mainly due to increased fourth 
stage error sensitivities. As an example the specific impulse sensitivity increased 
from 18.3 ft /sec (3 ) to 132.5 ft /sec. The coning rate effect also increased from 
12.3 ft/sec to 107 ft/sec. This increase is partly due to the fact that the errors 
have been propagated to 400,000 ft whereas in the polar case, they were depicted 
at fourth stage burnout. While most of the significant error sources showed some 
increase, the new fourth stage data and pitch profile has evidently enhanced the 
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effects of the fourth stage errors. The dispersion in flight path angle and 
inclination also increased from .262 deg and .258 deg respectively to .553 deg 
and .432 deg (la) due to the increased effect of coning rate errors. 

It should be noted that the deviations due to a variation in G (CMAL) 

ma 

of .6% (3a) was 8 ft and .07 ft/sec. The proper value to use was • 15% but an 
increase in the values by a factor of 25 will still be negligible. 



TABLE 2. NOMINAL REENTRY TRAJECTORY PARAMETERS 


Time (sec) 

0 

3 

8 

31 

44 

101 

183 

280 

305 

Coast Periods: 

First Stage Coast - 2.5 sec 
Second Stage Coast - -142.22 sec 

Third Stage Coast - 21.00 sec 

Launch Azimuth - 127*35 deg 
Payload Weight - 476 lbs 
Launch Position 

Latitude (geodetic) - 37*8479 deg 

Longitude - -75*4739 deg 
Altitude - 0 ft 

Fourth Stage Data: 


of Burn 

CG Location 
Station (inches) 

0 

42.7 

20 

41.2 

40 

39.2 

60 

37.1 

80 

33.8 

100 

29.4 


Pitch Rate (deg/sec) 
+.0265 
-3.00517 
- .81401 
-.475 

-.37069 
-.18193 
- .20677 

-.15 

-.2 


I (slug-ft 2 ) I (slug-ft 2 ) 

y x 


131.6 

8.62 

124.6 

8.16 

115.5 

7.50 

105.0 

6.57 

90.6 

5.42 

71.5 

4.01 
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TABLE 2. NOMINAL REENTRY TRAJECTORY PARAMETERS ( Coni' d) 


Spin Rate - 15 0 rpm 

Conditions at 400,000 ft 
Time - 506 . 87 sec 
Velocity - 26202 ft/sec 
Flight Path Angle - -3*065 deg 

Inclination - 50.425 deg 
Est. Time of Impact - 861.8 sec 
Est. Velocity at Impact - 26332 ft/sec 
Est. Latitude of Impact - 6.29 deg 
Est. Longitude of Impact - -43.6 deg 
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TABLE 3* SUMMARY OP THE STATISTICAL DISTRIBUTIONS, 
OF THE OUTPUT PARAMETERS 


OISPERSIGN OF THE SEMIMAUQ* AXIStFT. i 






11 ^ 04841 ^ 05 
flliii 8 T 49 E' 05 


2ND PERCENTILE SAMPLE 
5TH PERCENTILE SAMPLE 
PERC.ENJ I L 6 SAMPLE 



-3. 135257 5E 05 
-2.7359300E 05 
1.9867975E ft?. 


» 2 -. 56 < 



. f •*.?■>:! 


ECCENTRICITY DISPERSION 



mm piwir 


6 , 250709 2 D -04 
8 . 4873663 D -03 
- 3 . 8723297 E — 02 






2ND PERCENTILE SAMPLE 
5TH PERCENTILE SAMPLE 
95TH PERCENTI LE SAMPLE 


sample 


1.6123567E-02 
1.3393990E-02 
I .49149676-02 


l*951$%2Wrm c* : ? 

2.54851396-02 - " V . 





- . k ■ , ...■ '■■■ vftS* : 

■V •' -.*• ■■.-■.• ,.>• ';•• • ^ f- " 'l: ^ •- •-/ ' 

j-tv 'ini' ’• I • V ^7$2£'<?-£ : 4 ./ -.. •'• \ *•&$&& 1 


STANDARD DEVIATION 
SMALLEST SAMPLE 
2ND PERCENTILE SAMPLE 


llllifC/: 5TH PERCENTILE SAMPLE 

^ ' “* '■ 


4 .32764650— 01 
•1.1757078E 00 
•8. 059616 IE— 01 




wrnmmm 




- 6 .:. 6876984 E- 0 '^ 


LARGEST SAMPLE 


3 .41886046 00 


MEAN 


STANDARD DEVIATION 
SMALLEST SAMPLE 
2ND PERCENTILE SAMPLE 








5.64815450-01 
-4 .05663306 00 
-1.1476784E 00 









1.10696226 00 




iiSa 


LARGEST SAMPLE 


1.59038166 00 


gj '”'&i.S»jMENT 


, 'iil'P'ilTOii: 01 S PERS:TiMMT : 0'EGR6ES 1 


MMSiwwm 







||| 1 |*I 1 |^B PERCENTILE; SAMPLE |K|i| 
: >. : 95 TH percentile. SAMPLE ' ' * 


STANDARD DEVIATION 
SMALLEST SAMPLE 
2N D P ERCENT IL E S A MP LE 


7.57545660 00 
-8.4022601E 01 
- 1 m 269S 3 79 E . 0 1. 


WSSf 










r E 01 
1*215796 8E 01 

mmmzmmmttmMmm, Mm l 


LARGEST sample 


2.9031933E 01 


23 



DlSPERSfONOF arc LENGTH ALONG ORBIT (KM) 

MEAN » 2.0179964Q 00 


; Mfc a IM 

STANDARD DEVIATION 

3 

0*?64P 

2.243283 9D 

uu 

01 

SMALLEST SAMPLE 

3 

-6.3985839E 

01 

2ND PERCENTILE SAMPLE 

3 

-4.6811279E 

01 

5TH PERCENtHE SAMPLE 

• .. 3 

-3 .GR78784E 

01 

9 5TH PERCENTILE SAMPLE 

*' 

3.5184326E 

01 

98TH PERCENTILE SAMPtE 

• •-« 

4.1909301F 

01 

LARGEST SAMPLE 

3 

1.4571008E 

02 

RADIUS VECTOR DISPERSION 

(FT) 



MFAN 

3 

—5 .59 9781 30 

03 

STANDARD DEVIATION 

3 

4.5084073D 

04 

SMALLEST SAMPLE 

3 

-1.3468700E 

05 

2ND PERCENTILE SAMPLE 

3 

-1 .019345 OE 

05 

5TH PERCENTILE SAMPLE 

3 

-8.0582750E 

04 

95TH PERCENTILE SAMPLE 

; 3 

6.8894500E 

04 

98 TH PERCENT ILE SAMPLE 

3 

8.3900750E 

04 

LARGEST SAMPLE 


2.8475I50E 

05 

INERTIAL VELOCITY DISPERSIC 

m 

IF P.S). ; 


MEAN H- 

3 

— 1 .52903770 

01 

STANDARD DEVIATION 

3 

8 • 7046162D 

01 

SMALLEST SAMPLE 

3 

-4.6214233E 

02 

2ND PpRCEN T ILE SAMPLE 

3 

— 1 #8284058E 

02 

: 5TH PERCENTILE SAMPLE 

: 3 

-1.5982788E 


95TH PERCENTILE SAMPLE 

3 

I.2659253E 

02 

: 98TM PERCENTILE SAJMPLE 

m 

l .4581518E 

02 

LARGEST SAMPLE 

3 

2. 3220508E 

02 

AIRSPEED DISPERSION <FpS) 




MEAN - ’ 

. 3 

^l.SR47M#D 

01 

STANDARD DEVIATION 

3 

8 .77537320 

01 

SMALLEST SAMPLE 

3 

-4.7705200E 

02 

2ND PERCENTILE SAMPLE 

3 

-1.8989673E 

02 

, . 5TH PERCENTILE SAMPLE 

•3 

-1.61 16284E 

02 

95TN PERCENTILE SAMPLE 


1 .28 8064 OE 

02 

98TH PERCENTILE SAMPLE 

■ ■ 3 

1.4927026E 

02 ; : 

LARGEST SAMPLE 


2.3436499E 

02 

INERTIAL FLIGHT PATH ANGLE 

DISPERSION (DEGREES) 

MEAN 

3 

7.36498990-02 

STANDARD DEVIATION 

3 

5.5346475D- 

-01 

SMALLEST SAMPLE 

3 

-3.8355074E 

00 

2ND PERCENTILE SAMPLE 

3 

-1.0032183E 

00 

5TH PERCENTILE SAMPLE 

S 

-7.90 1788 IE-01 

95TH PERCENTILE SAMPLE 

3. 

9 .90 633 5 IE— 01 

98TH PlRCENTILi SAMPLE 

■ 

1 • 1529894E 

00 

LARGEST SAMPLE 

. 3 

1.7027636E 

00 
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ATMOSPHERIC fllGHT PATH ANGLE DISPERSION (DEGREES! 


5. 745668 OD-Ol 
3.9813019E 00 
1 .0409977E QQ 
8. 205503 5E-01 
1.Q285264E 00 


STANDARD DEVIATION 
SMALLEST SAMPLE 


5TH PERCENTILE SAMPLE 
95TH PERCENTILE SAMPLE 


LARGEST SAMPLE 


STANDARD DEVIATION 
SMALLEST SAMPLE 


7. 8277130E 01 
6.7718170E 01 


5TH PERCENTILE SAMPLE 
95TH PERCENTILE SAMPLE 


LARGEST SAMPLE 


PERIGEE OISPERSION IN Ml 


3.53495800 01 
1.4029141E 02 
8.S4275 2 0 E 01 
6.7454193E 01 
4.6378967E Ol 


STANDARD DEVIATION 
SMALLEST SAMPLE 




LARGEST SAMPLE 


PERIOD DISPERSION (SEC I 


5.38360840 01 
2.6143536E 02 


STANDARD DEVIATION 
SMALLEST SAMPLE 


1 .020596 9E 02 


: : : *iNf’PliCiN7fl^iSAWPL« 
9STH PERCENT HE SAMPLE 
98TH PERCENT 1 1 E S AMPLE 
LARGEST SAMPLE 


1.27IH69E 02 


IWGITDOE DISPERSION (DEGREES! 


STANDARD deviation 
SMALLEST SAMPLE 


2. 121906 3E-01 
1.60X0284E*01 


5TH PERCENTILE SAMPLE 
95th PERCENTILE SAMPLE 


LARGEST SAMPLE 
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TABLE 4. SIGNIFICANT ERROR SOURCES AND THEIR 
THREE SIGMA CONTRIBUTIONS AT REENTRY 


Error Source 

la Mag. 

Altitude 

(feet) 

Velocity 

(ft/sec) 

Fit. Path 
Angle 
(deg) 

ISP 

Specific 

.18% 

Impulse - First 

17084 

Stage 



MFRl 

Flowrate 

1.4% 

- First Stage 

-42371 

28.2 

-.296 

DKSG 

Torquer 

.35% 

Scale Factor 

-6434 

25.0 


DPBE 

3*57 mrad/ sec 

-17.4 



Pitch Rate Gyro Bias 

TMP1 1*67 mrad 11222 -113.0 

Pitch Thrust Misalignment - First Stage 

TMY1 1*67 mrad -8154 

Yaw Thrust Misalignment - First Stage 


CA01 1% -6801 

Drag coefficient - First Stage 


DRHO 

6.67% 

-44285 

22.2 

Density 




CDV1 

10% 

-4859 


Jet Vane Drag 



TIMl 

*078 sec 

11210 

-48.8 

First Timer 

Switch Uncertainty 


ISP2 

.094% 

4102 

15.2 

Specific Impulse - Second 

Stage 


MFR2 

1% 



Flow Rate - 

Second Stage 

13575 


TMP2 

1*67 mrad 

-19234 

49.9 


Pitch Thrust Misalignment - Second Stage 

TMY2 1*67 mrad 

Yaw Thrust Misalignment - Second Stage 

ISP 3 ,14% 33*8 

Specific Impulse - Third Stage 

MFR3 1*8% 87*1 

Flow Rate - Third Stage 


Inclination 

(deg) 


.216 



TABLE 4. (CONTINUED) 


Error Source la Mag. Altitude 

(feet) 

Velocity 
(ft/ sec) 

Fit. Path 
Angle 
(deg) 

Inclination 

(deg) 

TMP3 .557 mrad -32652 

Pitch Thrust Misalignment - Third Stage 

44.0 

-.336 


TMY3 .557 mrad 

Yaw Thrust Misalignment - Third Stage 



.254 

ISP4 .6% 

Specific Impulse - Fourth Stage 

132.5 



MFR4 1.8% 

Flow Rate - Fourth Stage 

15.08 



W4CP .03 rad/sec 11703 

Coning Rate - Pitch Axis 

29.3 

-.187 

1.01 

W4CY .03 rad/sec 96393 

Coning Rate - Yaw axis 

107.5 

1.38 

-.114 
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TABLE 5 


INDIVIDUAL THREE SIGMA ERROR SOURCE 
CONTRIBUTIONS AT SECOND STAGE IGNITION 
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TABLE 6. INDIVIDUAL THREE SIGMA ERROR SOURCE CONTRIBUTIONS AT 
THIRD STAGE IGNITION 
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TABLE 7. INDIVIDUAL THREE SIGMA ERROR SOURCE CONTRIBUTIONS AT 
FOURTH STAGE IGNITION 
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TABLE 10. DATA FROM CROSS TERM CARDS 
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KLO 5 XID b XLO 5 XtQ 5 XtQ b XLO b XLQ 

CUMULATIVE DISTRIBUTION OF DISPERSION OF THE SEMIfWOR AXIS (FT.? 



















































































































CUMULATIVE DISTRIBUTION OF ECCENTRICITY DISPERSION 
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CUMULATIVE DISTRIBUTION OF INCLINATION DISPERSION 4l>EGREEC 
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CUM'JLRTIVIE DISTRIBUTION OF LONG. OF ASCENDING NODE DISPEFTIuN CDEGREES? 
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CUMULATIVE DISTRIBUTION OF ARGUMENT OF PERIGEE DISPERSION (DEGREES? 
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CUMULATIVE DISTRISJTIuN OF DISPERSION OF ARC LENGTH ALONG PR=!TT kNM? 
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XIO 6 XIQ S XtO 5 XtO 4 XLQ * <U) XLO XU “ XtO 

CUMULATIVE DISTRIBUTION OF RADIUS VECTOR DISPERSION srf? 
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CUnULRTIVE DISTRIBUTION OF RIR3PEED DISPERSIO CFPS? 
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CUMULATIVE DIoTHlGUTla' OF Ir EhTIAL FLIGHT PATH AMGLE riFPERCION ^DEGREES} 
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curf LRTP'E DIJiRIBUTION OF PiFOSPHERi; FLIGHT rPTH PNGLE DIJPE^oICN ^TEGriEE 
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cum ILflTP'E DI3 TRT-3 JTJQi I OF APOGEE FISPERSION CNrD 
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PLIMULPTIVE riSTRieUTICM OF PERIGEE PISPERSION CNfl? 
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CUriULflTI'/E DIJTRIRUTIOi I OF PERIOD DI3PEK3ICM CSECP 
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nmULRTr'E DISTRIBUTION OF LONGITUDE DISPERSION -.'DEGREES? 
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CUMULATIVE DISTRIBUTION CF LATITUDE DISPERSION CDGREES? 
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CUMULATIVE DISTRIBUTION OF ALTITUDE DISPERSION CNM? 
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CUMULATIVE DICTRIBU riO, ' LF INERTIAL VERITY DICPERCION <Ff S? 



































ESCAPE MISSION ERROR ANALYSIS RESULTS 


The input data and the resulting nominal trajectory are listed in Table 11. 

In Table 12 are listed the means, standard deviations and the extreme values 
for each of the output variables. 

Table 13 lists the most significant error sources and their respective 
three sigma contributions to altitude, velocity, flight path angle and inclination. 

Tables 14, 15, 16, 17 and 18 are the individual three sigma contributions 
of each error source to altitude, velocity, flight path angle and inclination 
measured at each stage ignition and fifth stage burnout. 

Tables 19 and 20 lists the non-linear and cross-term data in this analysis 
in the same format as previously discussed. 

Figures 20-36 depict the cumulative distribution functions for the output 
variables. 

The results of the Monte Carlo analysis show that significant individual 
contributions are present in all of the stages with no signle source predominating, 
leading to a total velocity dispersion of 78 ft/sec. The inclination variation is 
less than in the four stage mission due to the decreased sensitivity caused by a 
large nominal velocity. 

In general, the five burn trajectory was less sensitive to errors in first 
and second stage than was the four stage mission. This effect compensated for the 
additional fifth stage errors. The large velocity bias effects present in the 
four stage mission was also decreased substantially. 
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TABLE 11. ESCAPE MISSION INPUT DATA AND NOMINAL 
TRAJECTORY PITCH PROFILE: 


Time (sec) 
0 

3. 

10 . 

31. 

42.2 

76.0 

83.42 

126.44 

144.83 

183-73 

Coast Periods 

First Stage Coast 
Second Stage Coast 
Third Stage Coast 
Fourth Stage Coast 

Launch Azimuth 


Pitch Rate (deg/sec) 
0 

-2.625 

- .813 

- .700 

- .461 

- .081 

- .432 
+ .400 

- .346 

- .569 


7.42 sec 

22.39 se c 

3.33 see 
9*5 sec 

90 deg 


Launch Position: 
Geodetic Latitude 
Longitude 
Altitude 


37.8479 deg 
-75.4739 deg 
0.0 feet 
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TABLE 1L ESCAPE MISSION INPUT DATA AND NOMINAL 
TRAJECTORY PITCH PROFILE: (Cont'd) 


Fifth-Stage Data: 


time- sec 

lbs Cons Wt remaining 

lbs Thrust 

Sta. c.g.-in 

Roll, l x 

Pitch ( y 

0.00 

194.85 

0.00 

15.36 

2.62 

12.57 

0.82 

179.35 

5905.0 

14.88 

2.56 

12.23 

1.67 

162.80 

608 O.O 

14.40 

2.46 

11.88 

2.53 

145.70 

6455.0 

13.92 

2.36 

11.54 

3.42 

128.23 

6210.0 

13.32 

2.24 

11.18 

4.33 

110.98 

6030.0 

12.72 

2.10 

10.86 

5.19 

93.93 

5975.0 

12.12 

1.97 

.10.62 

6.07 

76.85 

5960.0 

1 1 .40 

1.80 

10.36 

6.95 

59.94 

5940.0' 

10.80 

1.66 

10.12 

7.8] 

42.99 

5920.0 

9.96 

1.48 

9.89 

8.69 

26.00 

5600.0 

9.00 • 

1.30 

9.65 

9.55 

9.13 

50.0 

7.92 

1.10 

9.42. 

10.01 

. 0.00 

0.0 

7.56 

.98 

9.42 . 


Nominal Conditions at Fifth Stage Burnout 

Time : 240 sec 

Altitude: 679,627 ft 

Velocity: 39>199 ft/sec 

Flight Path Angle 11.77 deg 

O 

Semi-major Axis -60.673 x 10 ft 
Eccentricity 1.342 
Inclination 37*66 deg 
Inertial Longitude 293*59 de g 
Declination 37*317 deg 
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TABLE -11. STATISTICAL SUMMARY OF THE ESCAPE MISSION 
ERROR ANALY 




DISPERSION OF THE SEMIMAJOR AXIS (FT.) 


, V, 




STANDARD DEVIATION 
SMALLEST SAMPLE 


2ND PERCENTILE SAMPLE 
5TH PERCENTILE SAMPLE 
95TH PERCENTILE SAMPLE 
98TH PERCENTILE SAMPLE 
LARGEST SAMPLE 


-5.40986940 04 
1.81691380 06 
— 5.0291464E 06 


-3. 9395505E 06 
-2.9435135E 06 
2.7749790E 06 
3. 5809730E 06 
6.0420049E 06 


HIM IlilMMliM 

;, v ; * ' s ' : ''t-M'M 



INCLINATION DISPERSION (DEGREES) 




STANDARD DEVIATION 
SMALLEST SAMPLE 
2ND PERCENTILE SAMPLE 
5TH PERCENTILE SAMPLE 
95 TH PERCENTILE SAMPLE 
98TH PERCENTILE SAMPLE 
LARGEST SAMPLE 


1. 6885141D-02 
4. 5 528 176 D-0 2 
-1. 8545580E-01 
-7.0724010E-02 
-5.6 1962 13E-02 
8.7185383E-02 
1.0744667€“01 
2.33647 82E-01 


llllli 



OF ASCENDING NODE DISPERSION (DEGREES) 


LONG 


3. 6453081D“02 
5.04302210-01 


STANOARD DEVIATION 


SMALLEST SAMPLE a -3.2027664E 00 



2ND PERCENTILE SAMPLE 


STM PERCENTILE SAMPLE 


8.6 25450 1 E— 0 1 
1.0262108E 00 
1. 5043888E 00 


95 TH PERCENTILE SAMPLE 
98TH PERCENTILE SAMPLE 
LARGEST SAMPLE 



STANDARD DEVIATION = 7. 6414563D-0 1 

SMALLEST SAMPLE » -2.2351227E 00 

2ND PERCENTILE SAMPLE » -1.5841312E 00 

5TH PERCENTILE SAMPLE * -1.2836466E 00 

95TH PERCENTILE SAMPLE * 1.1474323E 00 

98TH PERCENTILE SAMPLE « 1.4977131E 00 

LARGEST SAMPLE * 2.7279234E 00 
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LARGEST SAMPLE = 2.5061621E 02 



STANDARD DEVIATION = 7.8945611D 01 

SMALLEST SAMPLE * -2.3756836E 02 


; 

sfiifi’i 8 II ’ . ; 


2ND PERCENTILE SAMPLE 
5TH PERCENTILE SAMPLE 
V51H PERCENTILE SAMPLE 
98TH PERCENTILE SAMPLE 
LARGEST SAMPLE 


-1.5806348E 02 
«t# 30378901 02 
i« 2856494E 02 
1.6292920E 02 
2.4808057E 02 



INERTIAL FLIGHT PATH ANGLE DISPERSION I DEGREES! 


7 #23907540-03 


MEAN 


3.7744043D-01 


STANDARD DEVIATION 


SMALLEST SAMPLE * -9. 6414053E-01 


2ND PERCENTILE SAMPLE * -7. 5141764E-01 

5TH PERCENTILE SAMPLE * -6# 7253959E-01 

mmmmmi 95TH PERCENTILE SAMPLE * 6.2411928E-01 

__ 98TH PERCENTILE SAMPLE * 7.4625742E-01 

LARGEST SAMPLE » 1.4544786E 00 
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ATMOSPHERIC FLIGHT 

PATH ANGLE 

DISPERSION (DEGREES) 

MEAN 


safe 

-7.442478GD-03 

STANDARD DEVIATION 

9 

3,89576680-01 

SMALLEST SAMPLE 


9 

-9.9504232E-01 

2ND PERCENTILE 

SAMPLE 

9 

-7, 74781 58E-01 

5TH PERCENTILE 

SAMPLE 

as 

-6,9403 32 7E-01 

95TH PERCENTILE 
98TH PERCENTILE 

sample 

9 

6, 4494646E-01 

SAMPLE 

m 

7.6998055E-01 

LARGEST SAMPLE 


s 

1.5011804E 00 

v-..:, ; j APOGEE DISPERSION 

INMI 


-1. 79393280 01 

STANDARD DEVIATION 

9 

5, 9938976D 02 

SMALLEST SAMPLE 


- 

-I.6566I94E 03 

2ND PERCENTILE 

SAMPLE 

9 

-1, 30633I3E 03 

5TH PERCENTILE 
95TH PERCENTILE 
98TH PERCENTILE 

SAMPLE 

SAMPLE 

SAMPLE 

m 

9 

-9,71977041 02 
9.2136986E 02 
l. L822742E 03 

LARGEST SAMPLE 



l. 9788460E 03 

disp6rsion 

CNMI 

9 

1*32285840-01 

STANDARD DEVIATION 

9 

6. 3033234D 00 

SMALLEST SAMPLE 


9 

-1, 6904144E 01 

2ND PERCENTILE 

SAMPLE 

9 , 

-1, 2525665E 01 

percentile 

95TH PERCENTILE 
98TH PERCENTILE 

SAMPLE 

SAMPLE 

SAMPLE 

m 

* 

-1, 07426 I5E 01 
1.0363922E 01 
1.37882381 01 

LARGEST SAMPLE 


9 

2, 2240722E 01 

PERIOD DISPERSION 

(SEC) 



MEAN ; 


9 

-4.18818470 01 

STANDARD DEVIATION 

9 

1.1244949D 03 

SMALLEST SAMPLE 


9 

-3. 1755149E 03 

2ND PERCENTILE 

SAMPLE 

9 

-2.4768352E 03 

5TH PERCENTILE 

SAMPLE 

9 

*1,84327881 03 

95TH PERCENTILE 

SAMPLE 

SB 

1.6972954E 03 

98 TH PERCENTILE 

SAMPLE 

m 

2,18278526 03 

LARGEST SAMPLE 


9 

3.6439482E 03 

L 

3N (DEG RE 

ES) 


MEAN 


9 

7.14427890-03 

STANDARD DEVIATION 

9 

2,80989730-02 

SMALLEST SAMPLE 


9 

-7, 8540802E-02 

2ND PERCENTILE 

SAMPLE 

9 

-4. 8784256E-02 

5TH PERCENTILE 
95TH PERCENTILE 

SAMPLE 

SAMPLE 

9 
9 : 

— 4.0257454E-02 
5. 46 20743 E-02 

98TH PERCENTILE 

SAMPLE 

9 

6,57129291^02 

LARGEST SAMPLE 


9 

1.0338879E-01 
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LATITUDE DISPERSION IDGREES) 


WmMmimW&inem ;-v; 


I *80036470-02 

STANDARD DEVIATION 


4*45715650-02 

SMALLEST SAMPLE 

=- 

—2* 1982241 E-Ol 

2ND PERCENTILE SAMPLE 

s 

—7* 1645737E-02 

5TH PERCENTILE SAMPLE 

* 

—5*78956601—02 

95TB PERCENTILE SAMPLE 

■ s 

0*9 10942 IE— 02 

98TH PERCENTILE SAMPLE 


1 • 02 19049E— 01 

LARGEST SAMPLE 

3S 

1* 5267134E-01 


.... 

ALTITUDE DISPERSION (NMI 




; 

MEAN 

•as 

4*24979510-01 



STANDARD DEVIATION 

:= 

5*02922930 00 



SMALLEST SAMPLE 

= 

-1* 5371993E 01 



2ND PERCENTILE SAMPLE 

— 

-1.0243593E 01 


' : 

5TH PERCENTILE SAMPLE 
95TH PERCENTILE SAMPLE 

as 

-7.6180256E 00 
8.4205933E 00 


" : . " . • . . V‘ ' 

. 98TH PERCENTILE SAMPLE 

■ : :• as ■ 

1* 1 12182 5E 01 



LARGEST SAMPLE 

* 

1.6259771E 01 
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TABLE 13. SIGNIFICANT ERROR SOURCES AND THEIR THREE SIGMA 
CONTRIBUTIONS AT FIFTH STAGE BURNOUT 


Error Source la Mag. 

Altitude 

(feet) 

Velocity 

(ft/sec) 

Fit. Path 
Angle 
(deg) 

Inclination 

(deg) 

SIW1 .83% 

First Stage Inert Weight 

8720 


.0557 


ISP1 .18% 

Specific Impulse - First 

15973 

Stage 

18.07 

.1059 


MFR1 1.4% 

Flow Rate - First Stage 

18304 

19.71 

.1056 


DKSG . 35% 

Torquer Scale Factor 

6156 

24.36 

.0826 


TMP1 1.67 mrad 

Pitch Thrust Misalignment 

14741 

- First Stage 

-65.43 


.0639 

TMY1 1.67 mrad 

Yaw Thrust Misalignment - 

6895 

First Stage 




CA01 1% 

Drag Coefficient 

6615 




DRHO 6.67% 

Density 

4218 


.2673 


FWN1 

Wind Effects 




.0179 

TIMl .078 sec 

Timer Step Uncertainty 

12343 

28.96 

.0739 


ISP2 .094% 6401 

Specific Impulse - Second Stage 




MFR2 1% 

Flow Rate - Second Stage 

7002 

16 . 63 

.0596 


DBY2 10% 

Yaw Dead Band - Second Stage 



.0136 

TMP2 1.67 mrad 

-51023 

54.89 

-.0349 



Pitch Thrust Misalignment - Second Stage 

TMY2 1.67 mrad .0554 

Yaw Thrust Misalignment - Second Stage 

ISP 3 .14% 27.17 .0687 

Specific Impulse - Third Stage 
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TABLE 13. (CONTINUED) 


Error Source 1 Mag. Altitude 

(feet) 

Velocity 

(ft/sec) 

Fit. Path 
Angle 
(deg) 

Inclination 

(deg) 

MFR3 1. 8% 

Flow Rate - Third Stage 

24 . 44 



TMP3 .557 ffi-rad 14860 

Pitch Thrust Misalignment - Third Stage 

41.95 

.1194 


TMY3 .557 mrad 

Yaw Thrust Misalignment - Third Stage 



-.0136 

ISP 4 .6% 

Specific Impulse - Fourth Stage 

14.99 



W4CP .03 rad/sec 

Pitch Coning Rate - Fourth Stage 

17.41 

-.0633 

.0467 

W4CY .03 rad/sec 

Yaw Coning Rate - Fourth Stage 

69.88 

.5886 


MFR5 1,8% 

Flow Rate - Fifth Stage 

-52.66 

-.1054 


ISP 5 . 15% 

Specific Impulse - Fifth Stage 

48.94 



W 5 CP .0463 rad/sec 

Pitch Coning Rate - Fifth Stage 

-19.77 


.0719 

W5CY .0463 rad/sec 

-25. 71 

.5726 



Yaw Coning Rate - Fifth Stage 



TABLE 14. INDIVIDUAL THREE SIGMA ERROR SOURCE 

CONTRIBUTIONS AT SECOND STAGE IGNITION 
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TABLE 15 . INDIVIDUAL THREE SIGMA ERROR SOURCE 

CONTRIBUTIONS AT THIRD STAGE IGNITION 
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TABLE 17. INDIVIDUAL THREE SIGMA ERROR SOURCE 

CONTRIBUTIONS AT FI PTH STAGE IGNITION 
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TABLE 18 . INDIVIDUAL THREE SIGMA ERROR SOURCE 
CONTRIBUTIONS AT FIFTH STAGE BURNOUT 
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TABLE 19. DATA FROM NON-LINEAR CARDS 
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TABLE 20. DATA FROM CROSS TERM CARDS 




in 

. 

m 


4 

m 

eg 


m 


, 

in 


in 

4 


eg 

eg 


«4 



o 

o 


o 

o 

O 


o 


O 


o o 


o 

O 


o 



UJ 

UJ 


UJ 

UJ 

UJ 


UJ 


UJ 


LU! UJ 


UJ 

UJ 


UJ 



s* 

O' 


m 

o 

eg 


in 

!00 


eg 

eg 


O' 

f- 


eg 



* 

ao 



s- 

O' 


in 


gs 


O' 

g> 


go 

f-4 


eg 




20 


eg 

o 

eg 


m 


gj 


e - 1 4 


m 

m 


in 



—4 

eg 


in 

in 

m 


O' 


S* 


ao. ^4 


g> 

in 


m 



o 

4 


4 

4) 

g> 


O' 


m 


O' 


O' 

eg 


m 



o 

ao 


O 

4 

eg 


m 


4 


m 

4 


4 

gs 


e- 



• 

• 


• 

• 

• 


• 


• 


9 

• 


9 

• 


9 



eg 

m 



4 

O' 


—4 


eg 


g> 

in 


H 

<-4 


-4 



1 

8 


i 


« 


1 


1 


8 

8 





1 



m 

eg 


eg 

4 

in 


4 

!m 


mi eg 


4- 

4 


4 



o 

o 


o 

° 

°j 

O 


o 


° 

° 


o 

O 


O 



UJ 

UJ 


LU 

UJ 

J 

UJ 

. 

LU 


LU 

UJ 

;uj 

UJ 


UJ 



CO 

H 


m 

a* 

4 


ao 


h- 


o 

4 


00 


g> 




-4 


tn 

m 

s- 


e- 


0 ' 


s 0 

m 


o 

ao 


s- 



->4 

4 


*■4 

s- 

00 


g> 


s* 


mi m 


h» 

e- 


s* 



4 

o 


o 

m 

N* 


s- 


eg 


o 

ao 


go 

co 


4 



g> 

m 


m 

m 

in 


ao 

n 0 


O' 

4 


m 

in 


o 



in 

4 


4) 

4 

#4 


O 




m 

eg 


m 

m 


m 



• 

• 


• 

9 

• 


• 


• 


* 

9 


• 

• 


9 



CM 

g> 


ao 

m 

4 


eg 


m 


s* 

g> 


g> 

ao 


o-4 



1 

i 


1 






l 


I. 

8 








m 

4 


m 

in 

in 


4 


m 


4 

m 


in 

in 


4 



O 

o 


o 

° 

o 


O 


o 


O 

o 


o 

O 


O 



UJ 

LU 


UJ 

UJ 

LU 


UJ 


LU 


UJ 

LU 


LU 

tu 


LU 



O' 

m 


m 

in 

00 


ao 


4 


m 

m 


00 

CO 


00 



o 

oo 


e- 

s- 

-4 


o 


4 


i-4 

o 


4 

eg 


4 



m 

in 


h- 

>0 

O' 


s* 


m 


m 

e- 


O 

ao 


h» 



O 

4 


m 

in 

m 


o 


ao 


m 

ao 


o 



m 



o 

eg 


eg 

g> 

4 


O' 


O' 


m 

m 


eg 

CO 


00 


! 

m 

m 


-4 

m 

m 


m 


o 


o 

S“ 


eg 

; m 


s* 



• 

• 


• 

» 

9 


• 


• 


• 

• 


9 

9 


9 


j 

h- 

m 


eg 

eg 

in 


O' 


00 


4 



m 

eg 


ao 


i 

1 

i 


1 

1 

i 


i 




g 




i 




m 


eg 

m 

mi 

i 

eg 

m 


4 


eg 


4 

4i 

eg 

4 



o 

K O CO 

o 

s- o m o in 

o in 

ON O 


0 4 O 4 

O 

m 


l 

O 

1 o 

l 

o « a 

I O 

l 

o 

1 

o 

l 

o 

8 O 

8 

o 

1 o 

1 

O 


UJ 


LU 

UJ 

LU 


UJ 

UJ 


UJ 


LU 


LU 

UJ 


LU 

LU 



o 

UJ O UJ 

o 

LU O 

LU O LU 

o uj o 

uj O 

UJ 

O LU 

o 

LU 

O LU 

O UJ 


o 

So® 

o 

m o 

40<0 

o m o 

ao O -4 

O 4 

a 

vO 

O ce 

o 

oo 


o 

oo^ 

o 

00 O 

o o m 

o o o 

r4 o m 

O ^4 O 

eg 

os- 

o o 


o 

o o eg 

o 

o o 

4 o in 

o o o 

e- o 

00 

O eg 

o 

eg 

O 4 

o 

in 

(M 

s- 

m eg in 

s- 

in s- 

4 0 0 

r- 

m 


s* o 

eg 

s- m 

h» 

eg 

O O' 

s- 

n 0 

• 

g> 

n o m 

g? s- o 

m 4 ao 

gj eg m 

h* 

»g 

00 

m m 

in 

00 

4 eg 

in 

O' 

> (/> 

9 

*4 

9 H 

• 

*4 • 

oo 

• o 

« 


• 

m 

e 

o 

9 «“4 

9 

>0 

# e- 

9 

4 

UJ X 


• 

eg * 

HI 

• 4 

4 4 • 


• m 

• 

eg 

9 

in ® »tn 

9 

4 • 

in 

9 

o o c 


-4 

r-4 


4 

4 

eg 


eg 


m 



O' ; 




m 

• LU 



* 


8 


l 


1 


8 


1 

» . 





8 

Q H* 

.4 


m 

-4 

*4 


m 

4 


4 

4 


■H 

-4 


4 

-4 



o 

mom 

o m o 

oos 

OvOO 

m o g3 

o m o 

e- 

OS- 

Q 

g> 

1 - V> 

» 

o 

8 O 

1 

O 1 

o 

1 u 

1 

o 

1 

o 

8 

O 

8 O 

f 

o 

1 o 

8 

o 


LU 


UJ 

LU 

UJ 


LU 

LU 


UJ 


UJ 


UJ 

LU 


JJ 

UJ 



o 

UJ O UJ 

O 

LU O 

UJ O LU 

O UJ o 

UJ O UJ 

O IJJ 

o 

UJ 

O UJ 

o 

UJ 

OC 

o 

SO® 

o o o 

eg O O' 

o^to 

gs o 


o m o 

eg 

O in 

o s- 

•4 O 

o 

O O ''0 

ONO 

ho® 

0 4-0 

g> o 

g^ 

o mo 


o m 

o 

o 

• 

o 

mom 

o 

ao o 

'f o o 

o r- o 

^ O N 

o o o 

in 

o o 

o 

in 

> O 

o 

so s» m 

o 

>o O 

S- 4 

© • o 

n n 

20 

0-4 0 

4 

s- m 

o 

g> 

UJ 

o 

4 g> m 

in k o 

O' g> 4 

in in © 

m 

? 

in 4 o 

hi in 4 

m 

s 

Q 

• 

00 

• «~4 

• 

s- • 

20 

• o*4 

• 


• 


• 

9 m 

e 


• 4 

9 

• 

i-4 

• 

•4 * 

eg 

• H 

9 

H • 

eg 

• 

4 

• in 

9 

eg 9 

-4 

9 

m • 

eg 

9 

o 


m 

m 


4 

O' 

eg 


m 


g> 


•—4 

«-4 


—4 

i4 


m 

x 


I 

1 


1 






l 


8 

8 






1- 




















v> 


g> 

g> 


m 

s- 

ao 


s* 


>0 


gj 

in 


s- 

s* 


s- 


(NJ 

O CM O 

non 

O fM O 

mom 

o m o 

mom 

o 

m o 

m o 

eg 

UJ 

> 

X 

> 

UJ ec UJ 

> a. 

ZUJI 

a. 

uj a£ uj 

£ 

UJ 

£ 

>- 

UJ ct UJ 

X UJ 

a 

uu 

& UJ 

CL 

X 

UJ 

Q 

8- 

-* X. fM 


m h- 

o * O' 

»- oo i- 

4 * ^ 

h* m i- 

m 


h in 

O 


4 

m 


•-4 

O' 

g> 


*4 


4 


00 

oo 


o 

s- 


m 

O 

1 

O' 

4 


eg 

eg 

m 


in 


O' 


4 

in 


o 

00 


4 



eg 

eg 


O' 

00 

i-4 


m 


m 


O' 

ao 


o> 

in 


m 

•>4 


e- 

o 


4 

eg 

m 


m 


HI 


O' 

go 


O' 

ao 


m 

UJ 

9 

fM 

IfMMN 

eg 4 - eg 

m eg o 

eg oo m 

o m O' 

m eg m 

oo mo 

m g» 

>! 

• > * 

UJ 

• Q, 

• a • 

UJ 

• > 

• > 

ft 

uj • a 

ft &. 9 

UJ 

• 

s 

flDj 

^ x g* 

O 4* ffl 

4 X ** 

O r* da 


o m ao 

in 

SET) 

o -• 

u 

a 

i 

i- I 

<* 

i a 

• 

8 

de 

8 

Ol 



» 

oe 

q 


at 

i 




93 


=6\ =4. oTo s . 4. e. s. 

XLQ 6 XLQ B XLQ B XLO B XLO B XLQ B XLQ B XLO 

CUMULATIVE DISTRIBUTION OF DISPERSION OF THE SEMIMAJOR AXIS (FT.? 
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CUMULATIVE DISTRIBUTION OF ECCENTRICITY DISPERSION 
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CUMULATIVE DISTRIBUTION OF INCLINATION DISPERSION CDEGREES) 
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CUMULATIVE DISTRIBUTION OF LONG. OF ASCENDING NODE DISPERSION < DEGREES? 
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CUMULATIVE DISTRIBUTION OF ARGUMENT OF PERIGEE DISPERSION (DEGREES) 
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: -5<3. O EQ. 1QQ. Iso: 200. 

CUMULATIVE DISTRIBUTION OF DISPERSION OF ARC LENGTH ALONG ORBIT CNF1? 
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CUNULATIVE DISTRIBUTION OF INERTIAL VELOCITY DISPERSION <FPS? 
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CUMULATIVE DISTRIBUTION OF AIRSPEED DISPERSION <FPS? 
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CUnULflTIVE DISTRIBUTION OF INERTIAL FLIGHT PATH ANGLE DISPERSION (DEGREES? 
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CUMULATIVE DISTRIBUTION OF ATMOSPHERIC FLIGHT PATH ANGLE DISPERSION (DEGREES? 
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CUMULATIVE DISTRIBUTION OF APOGEE DISPERSION CNM> 
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CUMULATIVE DISTRIBUTION OF PERIGEE DISPERSION CNM? 
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CUMULATIVE DISTRIBUTION OF PERIOD DISPERSION < SEC? 
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CUMULATIVE DISTRIBUTION OF LONGITUDE DISPERSION (DEGREES? 
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CUMULATIVE DISTRIBUTION OF LATITUDE DISPERSION CDGREES? 





109 


CUMULATIVE DISTRIBUTION OF ALTITUDE DISPERSION ^NM? 









































































CONCLUSIONS 


The results of the analysis of two missions are presented in this report. 
These results indicate that, although the overall accuracy is not radically chang 
the effects of any individual error source is strongly dependent on the particula 
mission trajectory which is used. 

Some general comments comparing the reentry mission, excape mission and the 
polar orbital mission which was analyzed in Phase I of this contract are in order 
The least sensitive trajectory is the excape mission since the large fifth stage 
burn reduces the effects of the lower stages, particularly in inclination, while 
burn itself is of such short duration that the trajectory errors are small. 

The most sensitive is the reentry trajectory, partly due to the fact that 
the errors are propagated from burnout to reentry but also due to an increased fo 
stage sensitivity. 
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